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Multiplicity fluctuations are studied both globaly (in terms of high-order moments) and locally (in terms of small 
phase-space intervals). The ratio of cumulant factorial to factorial moments of the charged-particle multiplicity 
distribution shows a quasi-oscillatory behaviour similar to that predicted by the NNLLA of perturbative QCD. 
However, an analysis of the sub-jet multiplicity distribution at perturbative scales shows that these oscillations 
cannot be related to the NNLLA prediction. We investigate how it is possible to reproduce the oscillations 
within the framework of Monte-Carlo models. Furthermore, local multiplicity fluctuations in angular phase- 
space intervals are compared with Monte-Carlo models and with first-order QCD predictions. While JETSET 
reproduces the experimental data very well, the predictions of the Double Leading Log Approximations and 
estimates obtained in Modified Leading Log Approximations deviate significantly from the data. 



1. INTRODUCTION 

The shape of a distribution can be analyzed 
in terms of its moments H. In the case of the 
multiplicity of particles produced in particle col- 
lisions at high energies, such an approach is pref- 
cred both theoretically (generating function) and 
experimentally (stability) . 

The normalised factorial moment of rank q of 
the multiplicity distribution P n is defined by: 



E~= g "("-l)--(rc-g+l)fn 



(1) 



It expresses the normalized phase-space integral 
over the g-particle density function. If particles 
are produced independently, the multiplicity dis- 
tribution is a Poissonian (for which the factorial 
moments of all ranks are equal to one). If the par- 
ticles are positively correlated, the distribution is 
broader than Poisson and the factorial moments 
are greater than one. 

The cumulant factorial moments are obtained 
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from the factorial moments by: 



K„ 



9-1 



(3-1)! 



*—' m\(q — m — 1) 

m=0 vy ; 



(2) 



K q expresses the normalized phase-space integral 
over the g-particle correlation function, i.e., that 
function which describes the genuine correlations 
between q particles. 

Here, we shall use these moments, for full phase 
space, as well as in ever smaller intervals of it, to 
study the interplay between hard and soft physics 
in particle production. In doing so, we use pertur- 
bative QCD (and its anomalous dimension) not as 
a model to be tested, but as a description of the 
perturbative background at the onset of unknown 
non-perturbative hadronization phenomena. 

2. THE GLOBAL MULTIPLICITY DIS- 
TRIBUTION 

2.1. High-order correlations 

Since \K q \ and F q increase rapidly with in- 
creasing q, it is useful to define g] the ratio 
H q = K q /F q , which has the same order of mag- 
nitude over a large range of q. One can view this 
ratio as a cumulant factorial moment normalised 




Moment Rank q 

Figure 1. Predictions at different order of pertur- 
bative QCD for the H q as a function of the rank 
3 [2]. 

to the factorial moment. It reflects the genuine q- 
particle correlation relative to the g-particle den- 
sity. 

The H q have been calculated for the gluon mul- 
tiplicity distribution at different orders of pertur- 
bative QCD (| (see Fig. 1): 

• For the Double Leading Logarithm Approxi- 
mation (DLLA), H q decreases to as q~ 2 . 

• For the Modified Leading Logarithm Approx- 
imation (MLLA), Hq decreases to a negative 
minimum at q — 5, and then rises to approach 
asymptotically. 

• For the Next to Leading Logarithm Approx- 
imation (NLLA), H q decreases to a positive 
minimum at q — 5 and goes to a positive con- 
stant value for large q. 

• For the Next to Next to Leading Logarithm 
Approximation (NNLLA), H q decreases to a 
negative first minimum for q = 5 and shows 
quasi-oscillation around for q > 5. 

Assuming the validity of the Local Parton 
Hadron Duality hypothesis (LPHD), such be- 
haviour is also expected for the charged-particle 
multiplicity distribution. 

The Hq of the charged-particle multiplicity dis- 
tribution have been measured by the SLD collab- 
oration H for e + e~ — ► hadrons at y/s = Mz- 
They indeed observe (Fig. 0) a negative minimum 
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Figure 2. H q versus q measured by the SLD col- 
laboration [3]. 



for q = 5 followed by quasi-oscillations about 
for higher q. This result has been interpreted as 
confirmation of the NNLLA prediction. 

Earlier, the H q moments were calculated ||,[| 
from charged particle multiplicity distributions 
of e + e~ experiments between 22 and 91 GeV. 
Similar behaviour with a negative first mini- 
mum followed by quasi-oscillation about is ob- 
tained. Furthermore, the same behaviour has 
been observed in hadron collisions between 20 and 
900 GeV |,§ and in hA an AA collisions § . 

While it is tempting to conclude that the ob- 
served behaviour is a confirmation of the NNLLA 
prediction, it must be realized that such a con- 
clusion rests strongly on the validity of LPHD. 
Although the LPHD hypothesis is accepted for 
single-particle distributions and mean multiplic- 
ity, difficulties have been detected when reso- 
nances are involved |8|,|9|. Further, heavy quark 
decays will cause a modification of the multiplic- 
ity distribution which is unlikely to be accounted 
by LPHD |o). For these reasons, we consider the 
assumption of LPHD as highly questionable. 

Therefore, we present measurements of H q not 
only of the charged-particle multiplicity distribu- 
tion but also of the sub-jet multiplicity distribu- 
tions. By considering sub-jets, i.e., jets obtained 
with small value of jet resolution parameter, we 
can, in some sense, follow the development of par- 



ton to hadron by varying the jet resolution pa- 
rameter (and hence the energy scale). The vari- 
ation of the mean sub-jet multiplicity has been 
found |l(| to agree well with the analytic pertur- 
bative QCD calculation in the perturbative region 
(energy scale greater than about 1 — 2 GeV) and 
with Monte Carlo in both the perturbative and 
non-perturbative region. These sub-jet multiplic- 
ity distributions should therefore enable us to ex- 
tend this analysis to the energy scales of partons 
and to that of perturbative QCD, avoiding (or at 
least lessening) dependence on the assumption of 
LPHD. 

2.2. The data 

The analysis is based on data collected by the 
L3 detector O] at the energy of Z, corresponding 
to approximately 1.5M hadronic Z decays. Events 
are selected using a procedure based on the track 
momenta measured in the Central Tracking De- 
tector including the Silicon Micro-vertex Detec- 
tor. About 1M events survive this selection. 

The charged-particle multiplicity distribution 
was corrected for the detector acceptance and 
efficiency using a matrix [ Jl2[ determined from 
Monte-Carlo events generated by JETSET £u| 
and passed through detector simulation. Ef- 
fects from selection procedure and initial state 
radiation are corrected from detector simulation. 
Charged particles originating from Kg and A de- 
cay are left in the sample. 

The same procedure was applied to correct the 
sub-jet multiplicity distributions. 

To evaluate the errors of the H q , we generated 
10000 multiplicity distributions from the exper- 
imental one, by allowing a Gaussian fluctuation 
of each P(n). The widths of the Gaussians were 
given by the statistical errors of the P(n). The 
statistical error on the H q is taken as the standard 
deviation of the 10000 values of H q . 

2.3. Charged-particle multiplicity analysis 

The detector-corrected charged-particle multi- 
plicity distribution is shown in Fig. ffl and com- 
pared to that of JETSET. Reasonably good 
agreement is found. 

To be less sensitive to the large statistical fluc- 
tuations in the high-multiplicity tail due to a fi- 
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Figure 3. The fully corrected multiplicity distri- 
bution compared to JETSET 7.4 PS. 

nite data sample, the largest multiplicities (N c h > 
50) are removed. The ratios H q are shown in 
Fig. [| for both the initial charged-particle multi- 
plicity distribution and the truncated one. These 
results are in qualitatively good agreement with 
those of the SLD collaboration. 

In both the truncated and the non-truncated 
cases, the H q have a negative first minimum at 
q = 5 and quasi-oscillations around for higher 
q, but the amplitudes are slightly larger in the 
truncated case and are found to increase for fur- 
ther truncation (not shown). 

The H q calculated from JETSET 7.4 PS agree 
well with the data both, for the entire and the 
truncated distributions, when the generated mul- 
tiplicities are truncated in the same way as the 
data. The H q calculated from HERWIG 5.9 @ 
do not agree with the data. They also show an 
oscillatory behaviour, but with a negative first 
minimum at q = 7 and all minima and maxima 
shifted by 2 units of q. 

2.4. Sub-jet multiplicity analysis 

Since the behaviour of the H q agrees qualita- 
tively with that predicted in NNLLA, it is nat- 
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Figure 4. H q obtained from the charged-particle 
multiplicity distribution. 

ural to interpret this agreement as a confirma- 
tion of the NNLLA prediction, and SLD g in- 
deed gave this interpretation. However, JET- 
SET agrees with the data and HERWIG shows 
qualitatively similar behaviour, in spite of parton- 
shower modellings not being at the NNLLA order 
of pQCD and, therefore, not expected to oscillate. 
That they do, together with doubt on the LPHD 
assumption, suggests that we view the NNLLA 
interpretation with scepticism. 

The sub-jet analysis assumes that the sub-jet 
multiplicity distribution is related to the parton 
multiplicity distribution at the energy scale cor- 
responding to the jet resolution. By choosing a 
scale where perturbative QCD should be appli- 
cable (> 1 GeV), we should be able to test the 
pQCD predictions for the behaviour of the H q 
without the LPHD assumption. We define sub- 
jets using the Durham algorithm |15| , minimizing 
the hadronization effects by associating particles 
with large angles and small transverse energies. 

It has been shown p^ , pl]Jlt| in previous papers 
(see for review JL7]) that sub-jet mean multiplic- 
ities agree well with the NLLA prediction in the 
perturbative region, i.e., for values of y cu t greater 
than 10~ 3 which corresponds to transverse ener- 
gies greater than ~ 1 GeV. 

Since this study starts with the charged- 
particle multiplicity distribution, we have also 
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Figure 5. The soft gluon production rate as a 
function of tjq is compared to the NLLA predic- 
tion. Furthermore at the top of the figure, is in- 
dicated the evolution of the behaviour of H q with 
the energy scales. 

constructed sub-jets using only charged tracks. 
The charged-particle and sub-jets multiplicities 
are, of course, identical for sufficiently small val- 
ues of y cut (y C ut < 1CT 8 ). We have verified that 
using both charged and neutral particles gives 
similar results. 

The average sub-jet multiplicities in 2-jet and 
3-jet events (M2 and M3) are obtained for a set 
of values of cut-off parameter yo: 10 -7 < j/o < 
2/i, where y\ = 0.01 is the cut-off parameter 
which resolves 2-jet and 3-jet events. M 2 {yo) and 
Ms(yo) correspond to the number of gluons radi- 
ated at different energy scales (at different places 
in the parton shower evolution) for 2-jet and 3-jet 
events, respectively. 



We have calculated the ratio 



A/3 -3 

A/2-2 



(2/o), which 



corresponds to the ratio of the amount of soft 
gluon production, when the primary quarks and 
the hard gluon have been taken out. This ra- 
tio, which has the advantage of being indepen- 
dent of Os, depends only on the resolution pa- 
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Figure 6. H q for sub-jet multiplicities, (a) for 
y cut = 6 x 10- 6 , (b) for y cut = 1.5 x RT 5 , (c) for 
y cut = 1 x 1CT 3 . 

rameters j/o an d y\. The NLLA prediction ]lS| 
with N f = 3 |l^ is seen (Fig. |) to be in good 
agreement with the data in the perturbative re- 
gion. This result supports the assumption that 
we can perform the sub-jet analysis with sub-jets 
obtained from charged particles. 

The sub-jet multiplicity distributions used for 
the calculation of the H q have been corrected us- 
ing the same procedure as used for the charged- 
particlc multiplicity distribution (i.e., correction 
matrix and correction for ISR, Kg and A decays). 
The resulting H q are shown in Fig. pi for three 
different values of y cu t . 

A first minimum at q — 4 followed by quasi- 
oscillation about 0, is seen forj/ cut = 6 x 10 -6 (en- 
ergy scale ~ 100 MeV) (Fig. ga). This behaviour 
is qualitatively similar to that for the charged- 
particle multiplicity distribution, as well as for 
sub-jet multiplicity distributions with smaller 
2/cutj but both the position of the first minimum 
and of the subsequent maxima and minima are 
shifted to lower values of q. 

As y cut is increased, the negative first mini- 
mum near 5 disappears. For q < 5 each H q alter- 
nates between positive and negative values. For 
higher q, the oscillations begin to disappear. For 
y cut = 1.5 x 10~ 5 (energy scale ~ 160 MeV) these 
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Figure 7. H q for the charged particle multiplicity 
distribution obtained from Monte Carlo genera- 
tors using different methods of parton generation 
and different models of fragmentation. 

oscillations have disappeared (Fig. ||b). 

In the perturbative region (energy scale > 
1 GeV), the oscillations and the negative first 
minimum near 5 have disappeared (Fig. oc). In- 
stead, H q alternates between positive and nega- 
tive values for each q (mind the scale!). 

Although good agreement with NLLA is found 
for the production rates of soft gluons seen in the 
perturbative region, none of the predictions made 
at different orders of pQCD for the H q is observed 
at LEP I energy. Even though pQCD describes 
collective processes (like the average sub-jet mul- 
tiplicity) it seems to be unable to describe the 
shape of these distributions. Therefore, the mini- 
mum and oscillatory behaviour of the H q seen for 
the charged-particle multiplicity distribution and 
for very low values of y cn t appears unrelated to 
the behaviour of the H q calculated in NNLLA. 

However JETSET PS reproduces the observed 
behaviour at all values of y C ut, and HERWIG 
shows a qualitatively similar behaviour. 

2.5. Monte Carlo study of the oscillatory 
behaviour of the H q 

Given the success of JETSET and HERWIG 
and the failure of the NNLLA prediction seen in 



the previous sub-section, we attempted to find 
that aspect of the Monte Carlo generators which 
is responsible for the agreement. 

We varied several options in JETSET and stud- 
ied their influence on the behaviour of the H q . 
First we used different models of parton genera- 
tion using in all cases the Lund string fragmenta- 
tion: 

• no angular ordering in the parton shower; 

• 0(a s ) and 0(a?) matrix elements; 

• matrix element production of only qq. 

Next we considered the possibility that the ob- 
served behaviour could come from the fragmen- 
tation model. We used JETSET 7.4 with the op- 
tions of parton generation given above but with 
independent instead of string fragmentation. We 
also used HERWIG 5.9 which uses cluster frag- 
mentation. 

In all cases the H q have a negative first mini- 
mum near q — 5 and quasi-oscillations for higher 
q. Some very different examples are shown in 
Fig. 0. This Monte-Carlo study shows that one 
can reproduce the oscillatory behaviour of H q 
without the need for NNLLA of pQCD. 

We further found that the oscillations do not 
exist on the parton level of JETSET at 91 GeV, 
even when reducing the cut-off Qq to 0.7 GeV. 
They do, however, set in at ^/s = 190 GeV. 

2.6. Conclusions from the global distribu- 
tion 

The ratio H q of cumulant factorial to facto- 
rial moments has been measured for both the 
charged-particle and sub-jet multiplicity distribu- 
tions. A perturbative QCD calculation at the 
next-to-next-to-leading-logarithm order predicts 
H q to have a negative first minimum at q ~ 5 and 
quasi-oscillation around for higher q, thus con- 
firming previous results. This behaviour is indeed 
observed for the charged-particle multiplicity dis- 
tribution. However, our analysis reveals that this 
behaviour appears only for very small y cu t, cor- 
responding to energy scales < 100 MeV, far away 
from the perturbative region. Furthermore, an in- 
vestigation performed on very different models of 
parton generation and fragmentation shows sim- 
ilar oscillatory behaviour for all cases. 

Contrary to jj), we conclude that the oscilla- 



tions of H q in the charged-particle multiplicity 
distribution are unrelated to the behaviour pre- 
dicted by the NNLLA perturbative QCD calcula- 
tions at the given energy. 

3. LOCAL FLUCTUATIONS 

3.1. Introduction 

Local multiplicity fluctuations have been stud- 
ied for many years in terms of a variety of phase- 
space variables |lj , but only recently has substan- 
tial progress been made to derive analytical QCD 
predictions for these observables |2(],£l],g2| . As- 
suming LPHD once again, these predictions for 
the parton level can be compared to experimental 
data. DELPHI J23J has shown that the analyti- 
cal predictions tend to underestimate the correla- 
tions between particles in small angular windows 
for a QCD dimensional scale of A ~ 0.1—0.2 GeV. 
However, a reasonable agreement is achieved for 
an effective A ~ 0.04 GeV, surprisingly smaller 
than theoretical QCD estimates p4j. 

In an investigation p5[ of local multiplicity 
fluctuations in terms of bunching parameters pt| , 
L3 concluded that the local fluctuations inside 
jets are of multifractal type. This is qualitatively 
consistent with the QCD predictions. Here, we 
extend this study and present a quantitative com- 
parison of the analytical first-order QCD predic- 
tions p0| , pl| , p2[ with the L3 data using normalized 
factorial moments of ranks q — 2 to 5 in angular 
phase-space intervals. 

3.2. Analytical predictions 

QCD is inherently intermittent and QCD pre- 
dictions [£0 21 2^] for normalized factorial mo- 
ments (NFMs) have the following scaling behav- 



^^(| 



(D-D q )(q-1) 



(3) 



where O is the half opening angle of a cone 
around the jet-axis, is the angular half-width 
of a ring around the jet-axis centered at 0o (see 
Fig. N), D is the underlying topological dimen- 
sion (D — 1 for single angle 6), and D q are the 
so-called Renyi dimensions. The analytical QCD 
expectations for D q are as follows p1,Sl],E2) : 



jet axis 




Figure 8. A schematic representation of the in- 
terval in the polar angle around the jet axis. 

1) In the fixed coupling regime, for moderately 
small angular bins, 



D q = 70 (0) 



9 + 1 



(4) 



where 7o(Q) = \f2Cxa s (Q)/iT is the anomalous 
QCD dimension calculated at Q ~ EQq, E = 
y/s/2, and gluon color factor Ca = N c = 3. 

2) In the running-coupling regime, for small 
bins, the Renyi dimensions become a function of 
the size of the angular ring (a s (Q) increases with 
decreasing 9). It is useful to introduce a new 
scaling variable pi] , 



ln(9 Q /e) 
ln(£0o/A) ' 



(5) 



where the maximum possible region (0 = 0o) 
corresponds to z = 0. 

Three approximations derived in DLLA will be 
compared: 

a) According to [glj , the D q have the form 



D q ~MQ)- 



9+1/. g 2 + l 
1 H r-^— -' 



(6) 



b) Another approximation has been suggested 
in m: 



D q ~2 l0 {Q)- 



q + 1 (\-^n~ 



9 



(7) 



c) In [£0J , a result has been obtained for the cu- 
mulant moments converging to factorial moments 



for high energies, 

q-w(q,z) 



D q ~2 l0 {Q)- 



z(q-l) ' 



w(q,z) = q\/l - z I 1 



m(l - z) 
V 



(8) 
(9) 



an estimate for D q has, furthermore, 



In [£1| 

been obtained in the Modified Leading Log Ap- 
proximation (MLLA). In this case, (ra) remains 
valid, but 7o(<2) is replaced by an effective 7q S (Q) 
depending on q: 
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where nf is the number of flavors. 

For our comparison of data and theoretical pre- 
dictions quoted above, we will use rif = 3 and 
A = 0.16GeV. 

The low value of m is chosen since even at 
high energies the production of higher flavors will 
rarely happen in the jet and its evolution is still 
dominated by the light flavors p9| . (The contri- 
bution of heavy quarks will be discussed below.) 
The value of A chosen is that found in tuning the 
JETSET 7.4 ME program |TJ on L3 data |T§ 
and in our most recent determination of a s (mz) 
@. 

For the angle Oo, we consider two possibilities: 
0o = 25° and 35°. The first value is chosen in 
order to compare our results with the DELPHI 
analysis |p3f , while the larger value of ©o allows 
a larger range of 6 to be studied. 

The effective coupling constant is evaluated at 
Q ~ Ee - For 6 = 25°, one obtains a s (£:e ) ^ 
0.144 according to the first-order QCD expression 
for a s (Q). This value leads to 7o(£6 ) - 0.525. 
For 6 = 35°, a s (EO ) - 0.135 and 7o(£6 ) - 
0.508. 

3.3. The data 

The data as well as the selection and correction 
procedures are the same as in Sect. 2.2, except 



that a more severe cut is applied on the track 
quality and K®, A decay products are excluded 
from the sample. Furthermore, Bose-Einstein 
correlations and Dalitz decay are removed from 
the data using correction from MC. 

The resolution of the L3 detector for a number 
of relevant variables was estimated in J28J. The 
resolution of polar angle defined with respect to 
the thrust axis is found to be approximately 0.01 
radians. For higher-order NFMs, the minimum 
angle is chosen according to the many-particle 
resolution studied in Q . 

The error bars on the results include contribu- 
tions from statistical errors on the raw quantities 
and both statistical and systematic errors on the 
correction factors. 

As a systematic error, we take half of the dif- 
ference between the correction factors determined 
using JETSET f| and those using HERWIG 5.9 
|[14| . In addition, we checked how a variation of 
the cuts affects the NFMs. The influence of such 
variations on the observed signal was found to be 
negligible. 

The sphericity axis is used to define the jet axis. 

3.4. Analysis 

3.4.1. Comparison with JETSET on ha- 
dronic and partonic levels 

It has been suggested in |22J to consider the 
ratio F q (z)/F q (0) in order to reduce hadroniza- 
tion effects on the actual behavior of the NFMs. 
In addition, this reduces a theoretical ambiguity 
for the evaluation of NFMs in full phase space 
(z = 0). In terms of F q (z)/F q (0), the power law 
(0) reads as 
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F q (z) 
F q (0) 



z(l~D q )(q-l)lli 



EQp 

A ' 



(11) 



The behavior of the ln(F q (z) / F q (0)) as a function 
of z is shown in Fig. 0. 

To increase statistics, we evaluated the NFMs 
in each sphericity hemisphere of an event and av- 
eraged the results, thus assuming that the local 
fluctuations in each hemisphere are independent. 

The open symbols show the predictions of the 
JETSET 7.4 PS model for hadronic (open circles) 
and partonic (open triangles) levels with the de- 
fault cut-off Q = 1.0 GeV. The Monte-Carlo 
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Figure 9. The behavior of F q (z)/F q (0) (q = 
2, . . . , 5) and JETSET 7.4 PS predictions, on the 
partonic and hadronic levels. Data are corrected 
for detector effects, initial-state photon radiation, 
Bose-Einstein correlations and Dalitz decays. 

prediction for hadrons gives a reasonable descrip- 
tion of the fluctuations. 

The partonic level of JETSET is close to the 
data, both for Q = 1.0 GeV and Q = 0.7 GeV, 
though a noticeable difference between the slopes 
for the data and the partonic level of JETSET is 
present. The similarity of partonic and hadronic 
level Monte-Carlo predictions is used as a mea- 
sure of the degree of validity of LPHD. 

The contributions of heavy flavors (c and b 
quarks) were estimated by rejection of these fla- 
vors on generator-level of JETSET. The effect 
was found to be negligible. 

3.4.2. Comparison of the analytical predic- 
tions 

The comparison of the analytical QCD predic- 
tion s (ffl -JlQl) to the corrected data is shown in 
Fig. Hgfor A = 0.16 GeV and 6 = 25°. For sec- 
ond order, the running-a s predictions lead to the 
saturation effects observed in the data, but signif- 
icantly underestimate the observed signal. Pre- 
dictions for the higher moments are too low for 
low values of z, but tend to overestimate the data 
at larger z. The fixed coupling regime (thin solid 
lines) approximates the running coupling regime 
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Figure 10. The analytical QCD predictions for 
A = 0.16 GeV: a s = const (4); DLLA (a) (6); 
DLLA (b) (7); DLLA (c) (8); MLLA (10). 



for small z, but does not exhibit the saturation 
effect seen in the data. The DLLA approximation 
(||) and (m) differ significantly at large z. The pre- 
dictions for cumulants (8) show a stronger satu- 
ration effect. The MLLA predictions do not differ 
significantly from the DLLA result (H). 

At 0o = 35° (not shown) both data and predic- 
tions show a stronger increase than for © = 25°. 
This indicates that fluctuations are larger for 
phase-space regions containing a larger contribu- 
tion from hard-gluon radiation. 

A better agreement of the QCD predictions 
with the data at low z-values can be achieved 
by decreasing the value of A. A similar observa- 
tion has been made by DELPHI |23|| . As an ex- 
ample, Fig. [Ly shows the case of A = 0.04 GeV. 
Such an effective value makes the coupling con- 
stant smaller and this can expand the range of 
reliability of the perturbative QCD calculations 
(for A = 0.04 GeV, a s (EQ ) ~ 0.112, j (EQ ) ~ 
0.46). However, this leads to a large disagreement 
between the QCD predictions and the data for 
z > 0.3, where contributions from high-order per- 
turbative QCD and hadronization are expected 
to be stronger. We have varied A in the range of 
0.04 — 0.25 GeV and found that there is no value 
of A in this range which provides agreement for 
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Figure 11. Same as Fig. 10, but for A = 0.04GeV. 

all orders of NFMs. Furthermore, increasing ri{ 
to 4 and 5 reduces the second moment to repro- 
duce the data at large z, but high-order moments 
still continue to overshoot. 

3.5. Conclusions from local fluctuations 

The first-order predictions of the DLLA of 
perturbative QCD are shown to be in disagree- 
ment with the local fluctuations as observed for 
hadronic Z decay. This conclusion is valid for 
standard values of A (A = 0.16 GeV) as well as 
for small values (A = 0.04 GeV). In the latter 
case, a reasonable estimate for z < 0.3 can be 
reached, consistent with the DELPHI conclusion 
p3|. However, our analysis shows that, in this 
case, the theoretical NFMs strongly overestimate 
the data for relatively large z (small 0), where 
contribution from high-order perturbative QCD 
is stronger. 

The prediction of MLLA shows very similar 
a result as DLLA. Note, however, that a full 
MLLA correction for the z-dependence of D q has 
not been obtained, so far. The MLLA predic- 
tion quoted above only modifies 70, but the z- 
dependence of NFMs is parameterized by DLLA, 
which is only asymptotically correct. One large 
effect in the MLLA approximation is the angular 
recoil effect, which is important for small z. This 
effect can change the value of NFMs at z — and, 
hence, the absolute normalization of the NFMs 



p9| . Note that a recent study J3(| 0I energy- 
conservation in triple-parton vertices shows that 
the energy-conservation constraint is indeed size- 
able and leads to a stronger saturation effect. 

Another likely reason for their disagreement 
with the experimental data is the asymptotic 
character of the QCD predictions, corresponding 
to an infinite number of partons in an event. Fur- 
thermore, the failure of the predictions can lie 
with the LPHD hypothesis, which is used to jus- 
tify comparison of the predictions of perturbative 
QCD to hadronic data. 
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